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In the scalar-tensor theory of gravitation it seems nontrivial to establish if solutions of the cosmological 
equations in the presence of a cosmological constant behave as attractors independently of the initial values. We 
develop a general formulation in terms of two-dimensional phase space. We show that there are two kinds of 
fixed points, one of which is an attractor depending on the coupling constant and equation of state. In the case 
with a power-law potential in the lordan frame, we also find new type of inflation caused by the coupling to the 
matter fluid. 



I. INTRODUCTION 

Einstein's General Relativity has proved to be the simplest 
theory for successful understanding of a number of experi- 
ments and observations. Still, on the other hand, there seem 
to be growing indications that a yet-to-be-discovered scalar 
field might play fundamental roles in cosmology. A list of 
the possible sources include the one expected from higher- 
dimensional theories of gravitation such as superstring/M the- 
ory and a scalar field (a volume modulus) which couples to 
4-dimensional gravity after compactification, also those in the 
context of the brane world scenario. We may reasonably ex- 
pect that realistic consequences of these hypotheses can be im- 
plemented in terms of the scalar-tensor theory proposed first 
by Jordan J2], developed later by Brans and Dicke 0]. 

One of the recent focuses of the scalar field is aimed partic- 
ularly at the origin of the dark energy which appears to be re- 
quired from the observed acceleration of the universe, with the 
renewed interest of today's version of the cosmological con- 
stant problem, culminating to the twin questions; fine-tuning 
problem and the coincidence problem. We argued 0, 0, [H [g] 
that the scalar-tensor theory is precisely that causes the behav- 
ior, A ~ t" 2 realized numerically by 10~ 120 ~ (10 60 )" 2 0], 
where 10 60 is today's age of the universe in units of the Planck 
time ~ 10 _43 s, hence preparing another simple implementa- 
tion of the scenario of the decaying cosmological constant as 
discussed in 0] and 0]. 

Before reaching conclusions to be compared with observa- 
tion, however, we must go through certain complications in- 
cluding details on the choice of the conformal frames among 
other things. Also many consequences derive from leaving 
the Brans-Dicke model 0], as required ultimately by a single 
technical aspect on the attractor nature of the cosmological so- 
lutions. This motivated us to develop a general framework of 
studying dynamics of the system including the scalar field tak- 
ing the unique roles of the conformal transformation properly 
into account. 
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After a brief introduction of the action of the scalar-tensor 
theory in both conformal frames, the Jordan and the Einstein 
frames in Section II, we enter Section III to develop a for- 
mulation in the Einstein frame in which we may trace how 
the cosmological solutions evolve in two-dimensional phase 
space. We assume the presence of the exponential potential 
of the scalar field corresponding to the simple cosmological 
constant in the Jordan frame. It is crucially important to make 
a right choice of the new extended time-coordinate other than 
the conventional cosmic time. Also of central importance is 
to deal with self-autonomous systems. 

In Subsection III-A we discuss the fixed points in phase 
space. We show that there are two different sets of fixed 
points: one (FP1) is for the well-known universe of scalar- 
field dominance, and the other (FP2) represents a new type of 
the universe in which the matter fluid energy is scaled to the 
potential of the scalar field because of the coupling with the 
matter fluid. We find that the universe of the latter type al- 
ways expands in the same way as the radiation-dominant uni- 
verse in the Einstein frame. In Subsection III-B, we present 
the stability analysis and the behavior of the attractor solu- 
tions. Flows of trajectories in phase space are illustrated for 
the examples taken from the radiation-dominant universe. Ba- 
sically the same analysis will be repeated in Section IV now 
in the Jordan frame. 

In Section V we generalize the argument to the power-law 
potential with a monomial of the scalar field multiplied by the 
cosmological constant chosen in the preceding section in the 
Jordan frame. We find a new type of inflation at the fixed 
point FP2 even for the potential which is too steep to cause 
inflation without the coupling with the matter fluid. Section 
VI is devoted to the concluding remarks. 

In Appendix A, we discuss the effect of the curvature term 
not included in the preceding sections. In the subsequent three 
appendices, we add related discussions on the accelerating 
universe. Appendix B reveals the presence of a complication 
in the numerical analysis in the Jordan frame, while in Ap- 
pendix C we discuss the attractor nature of the scale-invariant 
model as an alternative to the Brans-Dicke model. The final 
Appendix D will be devoted to offering another simplified ap- 
proach to the power-law potential. 
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II. SCALAR-TENSOR THEORY: JORDAN FRAME VS 
EINSTEIN FRAME 

We discuss cosmology in the scalar-tensor theory of grav- 
itation. We assume the presence of a potential V(<p) in the 
Jordan frame. The action is 



where 



S — Sg + S m J 



where 



S„ 



5,1 



d i x^/—g 



d 4 x^gL m (g, cf>) , 



(2.1) 

(2.2) 
(2.3) 



with e = ±1 11011 . Note that the matter action S m is assumed 
to have no scalar field <j> according to the Brans-Dicke model 
J2l, in which Weak Equivalence Principle (WEP) is intended 
to be respected. The scalar-tensor theory defined in this way 
is equivalent to the traditional Jordan-Brans-Dicke theory with 
an added potential U((p) = V(<f>), often expressed as 



S a = 



-9 



^i?--(V^) 2 



U(<p) 



(2.4) 



with the Brans-Dicke constant uj = e/(4£) and ip = (£/2)(/> 2 . 

Although we may further extend this type of the scalar- 
tensor theory with an arbitrary function of the scalar field cf> 
multiplied with R, we confine ourselves to the original simple 
if, because it features global scale invariance except generally 
for the V{4>) term. 

The choice e = —1 in ( 12. 2t is closely related to string 
theory. The D-dimensional action for the zero-modes in the 
closed string sector is given ifTHl by 



d U Xy 



-2* 



i?G?)+4(v$r 



(2.5) 



which, re-expressed according to our own sign convention as 
in (1.30) of JH], corresponds to the first two terms in j2.2j 
with e = — 1 and £ = 1/4, hence uj = —1, by introducing 

= 2e-* 

We can alwa ys m ove to the Einstein frame by a conformal 
transformation Hi [Tl [Tl 



where 



with 



— O 2 

g^v ' g*i±v ^ g^v 

Q 2 = ^ 2 = exp(2Ccr) , 

c 2 ^ (6 + e r 1 ) _1 = (6 + 4 W )- 1 



(2.6) 



(2.7) 



(2.8) 



which defines a canonical scalar field a in the Einstein frame; 



S 




1 



2 2 

d 4 Xy/-g t L m (ip :t ,g*, a) , 



(V,a) 2 -K(<7) 



(2.9) 



K(ct) = exp (-4Cct) ^(0) 



(2.10) 



We mark the quantities in the Einstein frame with the subscript 
*, while those in the Jordan frame are left unmarked, unless 
otherwise indicated. This is in accordance with the notation 
used in yfl. 

We point out that e = — 1, apparently indicating a ghost 
nature of the non-diagonalized field <j>, is a real difficulty only 
if ( 2 turns out to be negative implying a negative energy for 
the diagonalized field a. We always assume the condition 



C 2 > o. 



(2.11) 



This can be obeyed even if e = —1 if £ > 1/6. Imposing £ > 
0, which we assume throughout this paper, due to the required 
positivity of the energy of tensor gravity, we find that e = 1 
allows any £ but with ( 2 < 1/6, while e = — 1 constrains 
£_> 1/6 and ( 2 > 1/6, as displayed graphically in Fig. 1 of 

The parameters in ( 12.5b gives £ 2 = (D — 2)/4 which is 
1/2 for D = 4. If e = -1 and £ = 1/6, we find ( 2 -> 
oo, implying no kinetic term in the Einstein frame, hence no 
degree of freedom. We do not consider this choice any further. 

Because we assume that no <f> field enters L m in the Jordan 
frame, we find that the energy-momentum of the matter fluid 
is conserved in the Jordan frame; 



VTt = , 



(2.12) 



for which WEP is respected. The energy-momentum tensor in 
the Einstein frame is obtained by 



7^=exp(-4Ca)T^, 
which is no longer conserved; 



(2.13) 



(2.14) 



where T* = Tf p . Note that the universal free-fall (UFF) is 
still maintained, as an expression of WEP. 



III. COSMOLOGY WITH A COSMOLOGICAL 
CONSTANT : ANALYSIS IN THE EINSTEIN FRAME 

We discuss cosmology in the scalar-tensor theory with a 
cosmological constant V — Vq. In this section, we discuss it 
in the Einstein frame, though the analysis in the Jordan frame 
is given in the next section. 



A. The Basic Equations and the Fixed Points 

The metric of isotropic and homogeneous universe is given 
by the FRW form: 



dsi = —dt^ + azds 



3 ! 



(3.1) 



3 



where ds 2 is the metric of maximally symmetric three- 
dimensional space with the curvature constant k = or ±1. 
The basic equations in the Einstein frame are 



Hi 



1 /l 



3 V2 
<9V* 



I 



a + 3P»<r + V 1 = C(p* - 3P*), 
do 



(3.2) 
(3.3) 



where iJ» = d*/a, P*, and are the Hubble expansion pa- 
rameter, the pressure, and the energy density in the Einstein 
frame, respectively. The dot implies djdt* throughout in the 
Einstein frame. Eq. J2. 14-b in the Einstein frame is then re- 
expressed as 



3P»(P* + /?*) = — (<j(p* — 3P*) . 



(3.4) 



Assuming the equation of state P* = (7 — we further 
re-express Eqs. (I3.2b -( I3.4I ) into 



Focusing on k = 0, we differentiate Eq. d3.91 l with respect 
to r„, to obtain 



2-7/2 , 7 



37, 



G a + - + 2Qo' J£» -"Kl . (3.12) 



Here we have used Eqs. (13.91 ) and d3.10t as well as the equa- 
tion 



(p./V.)' - -3 7 (M* - C^')(P*/K) 



(3.13) 



which is obtained from Eq. (13. lit and the definition (12.101 ) of 

V*. 

In the same way we put Eq. (13.101 1 into 



V* 



(3.5) 



i 2 3 V2 

a + 3P»a-4CV; = C(4-3 7 )p„, (3.6) 
+ 3jH*p* = -C(4 - 37)07?*. (3.7) 

We now introduce a new dimension-free time coordinate r» 



by 



<ir* = 2\/Vi t dt 



(3.8) 



We further introduce IK, = a'^/a* , where the prime is for a 
differentiation with respect to r*. We then put Eqs. ( 13.5b - (13.7b 
into the new form 



2 V K 



3?0' - C 



(4-37)p. 



4K 



3 7 JC»p* = -C(4 - 3 7 )ctV. 



(3.9) 

= 0, (3.10) 
(3.11) 



a" = -3JC*a' + (2a' 2 + l) + 3C(4 - 3 7 )J{ 2 , (3.14) 



where Eq. d3.9t has been used to obtain the last term on the 
right-hand side. A set of equations ( 13.121 ) and ( 13.141 ) gives a 
self-autonomous system. In fact, by introducing the variables 
x and y defined by x = a' and y — (~~ ll X t , we derive 



x' = ^ [2 7 x 2 -4xy + 4C 2 (4-37)?/ 2 + 7] (3.15) 



y' = — [-2(2 - 7 )x 2 + 16C 2 ^ - 12C 2 7 y 2 + 7] 



(3.16) 



By choosing x' = y' = 0, we find four fixed-points in this 
system; 



J 



FP1± : (x F ,y F ) = (x[ ±} ,y^) = ± f 
FP2± : (x F ,y F ) = (x { ±\y { ± ] ) = ± ( 



2^ 



V^8^ 2Cy3 3 8?7 

V7 

V2(2- 7 -2(4-3 7 )C 2 )' ^2(2- 7 -2(4-3 7 )C 2 ) / 

I 



V7 



(3.17) 
(3.18) 



The fixed points FP1± exist if C 2 < 3/8, while the fixed 
points FP2 ± exist if 7 > 4/3 or if 7 < 4/3 with C 2 < 
(2 - 7)/[2(4 - 37)] (equivalently,7 > 2(4C 2 -1)/(6C 2 -1)). 
For £ = 1/2, two types of fixed points coincide to each other. 
In Fig. Q] we show in which portion of the £ 2 -7 plane we find 



the fixed points. 

From d3.9t with k = 0, we obtain 



El 
V* 



2(6CV -x 2 )-l. 



which is constant at the fixed points. 



(3.19) 



4 



A/2 



1/4 3/8 



FIG. 1: The fixed points FP1± and FP2± exist in the horizontally- 
and vertically-shaded regions, respectively. 



At the fixed point FP1±, we find immediately 
(p*/K)fpi = 0, 
while at the fixed point FP2 + , we obtain 



2(4C 2 - 1) 



FP2 2 - 7 - 2(4 - 3 7 )C 2 ' 
This result is consistent with Eq. ( 13.13b , i.e. 



-37C(y - x) 



V* 



(3.20) 



(3.21) 



(3.22) 



the right-hand side of which vanishes at the fixed point FP2 + , 
showing that the ratio of the energy density to the potential is 
constant. 

Next we discuss the scale factor and scalar field at the fixed 
points. We have 

a' = x F , (In a*)' = (y F , (3.23) 

for fixed points. 

ForFPl + ,we have (xf,Vf) — giving 



2C 



a* = a*o exp 



(3.24) 



By inverting ( 13.8b and substituting from ( 13.23b we obtain 

4C 2 



= i*o ex P 

Substituting this back into (13.25b and (13.26b then yields 



1 / t, 
a = ttt In — 



Co 



a* — do 



1/(8C 2 



(3.26) 

(3.27) 
(3.28) 



In order to fix o"o> we g° back to the original equations of 
motion, i.e, setting = finding 



H* = -( -& 2 + K 
* 3 I 2 



(3.29) 



which gives the value of the scalar field Co at i* = t*o as 

3-8C 2 



exp (-4£(r ) 



6K 4 t 2 v 



(3.30) 



Eq. (13.28b shows that the solution with £ < l/(2\/2) gives a 
power-law inflation. 

For the fixed point FP2 + , we have (xf,Vf) = 
(4 +) , yi +) ), replacing (13371 and dl28b by 



a = 



V2(2- 7 -2(4-3 7 )C 2 ) 
CV7 



a* = a*o exp 



r, +cr , (3.31) 
, (3.32) 



V2(2- 7 -2(4-3 7 )C 2 ) 
respectively. The cosmic time in the Einstein frame is 

2CV7 



£* = t*o exp 
Hence we find 



V2(2- 7 -2(4-3 7 )C 2 ) 



<r = llnf^ 

2C V*o 



Co 



1/2 



(3.33) 

(3.34) 
(3.35) 



1 /2 

It is important to notice that a* ~ tj follows also for dust- 
dominance. In fact this behavior is true for any equation of 
state. For the remedy of this unfavorable result, the reader is 
advised to see Section 4.4.3 of JH or Section 3.4 of 0]. 
To fix ctq, we use the original Friedmann equation: 



#2 = V + K / £ 

* 6 3 V* 



(3.36) 



v ' 7 where p*/V* is a constant given by Eq. ( 13.21I ). We then find 



» < At \ (4C 2 -l)[2-7-2(4-3 7 )C 2 ] n _ 

exp -44 (T = 1ft 7472 T/ • (3-37) 

I67C i*o^o 

Note that spacetime in this range is static in the Jordan frame, 
as will be shown in (I4.22b . 



B. Stability Analysis and the Attractors 

Next we analyze stability of the fixed points FP1± and 
FP2± in the self-autonomous system (13. 15b and ( 13.16b . 
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1. perturbation analysis 

The simplest way is to apply a perturbation analysis. We 
perturb the variable (x, y) around the fixed point (xp, yp) as 



(3.38) 







( Xp - 


\-6x\ 






V vf - 


V5y j 



Inserting Eq. ( 13.381 ) into the basic equations ( 13.15b and ( 13. 161 ). 
we find a set of the linear differential equations: 



For FP2 ± , we have 

(TrA) FP2± = -y { 2 ±} (3.48) 

(det^) FP2± = -37C 2 (3-8C 2 )(4C 2 -1)(^ ±} ) 2 (3.49) 

The stability condition for the expanding universe (the fixed 
point FP2 + ) gives 

1 



c 2 > 



(3.50) 



Sx \ _ / A xx A xy \ f 5x 



In this way we find that FP1 + and FP2 + are the attractors for 



5y J \ A yx A yy ) \5y 
where the components of the matrix A is given by 

' A xx = 3C("fx F ~ y F ) 

3C(-zf + 2C 2 (4-3 7 )2, f ) 



(3.39) 



4/3 



.2- 



(3.40) 



xf + 2(y F 



>-yx ~ 2 £ 

, A yy = C( 2x F ~ ZiDf) 



Assuming Sx, Sy oc e WT *, we find the equation for the eigen- 
value u> as 



' No attracto 



/ 



TiA lu + det A = , 



(3.41) 



where 



1/4 



FIG. 2: The attractor fixed-points (FP1+ and FP2+) in the same pa- 
rameter space of £ 2 and 7, also the same shading pattern as in Fig. 



□ 



+ {3 7 (4C 2 - 1) - (3 - SC 2 )}(x F - y F )] , (3.42) 

dGtA = 2(4^- 1) [(2 ~ 7 ~ 2(4 ~ 37)C2) {XF ~ 

- 2jC 2 (3-8(; 2 )(xF-yF) 2 ] . (3.43) 



C 2 < 1/4 and for C 2 > 1/4 (also with ( 2 < (2 - 7)/[2(4 - 
3 7 )] for the existence of FP2 + ), respectively. For C, 2 = 1/4, 
two types of fixed points merge with each other, sharing the 
same behaviors, but like a saddle point rather than an attractor. 
Fig. [2] shows in which portion of the ^-7 plane we have the 
attractor fixed points. 



In this expression, the first terms in (13.42l i and ( 13.431 ) vanish 
for FP1±, while the second terms disappear for FP2±. The 
fixed point is stable in the following two cases: 

• Eq. ( 13. 4U has two negative real roots. 



• Eq. ( 13. 4U has two complex conjugate roots with a neg- 
ative real part. 

The condition is 

TrA < 0, and dot A > . (3.44) 
For FP1±, we find 

(TrA) FP1± - [4(3 7 + 2)C 2 -3( 7 + l)]^ ±) (3.45) 



(detA) F pi 



3 7 C 2 (3-8C 2 )(4C 2 -l)( 2 /i ±) ) 2 (3.46) 



For the expanding universe (FP1 + ), which we are interested 
in, the above condition d3.44t gives 



c 2 < 



(3.47) 



2. phase space analysis 

We also study stability by use of a phase-s pac e analysis of 
the dynamical system with d345b and d3~T6l l 111 [Trjl [ni [H . 
We may discuss global stability rather than local one in the 
perturbative approach, as we will show shortly. There is a lim- 
itation, however, because dependence on the values of 7, £ 2 is 
not as simple as shown in ( 13.451 ) and in the subsequent equa- 
tions. We must develop the computation for each of these 
parameters separately, though without any difficulty in prin- 
ciple. For this reason, we show the following examples of 
radiation-dominance (7 = 4/3) illustrating generic features 
shared commonly by this type of analyses with any values of 
the parameters. 

For the sake of convenience we start with reproducing 
(13451 ) and (13461 ) for 7 = 4/3; 

x = ( (2x 2 - 3xy + 1) , (3.51) 
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which determine how the point (x, y) representing the solution 
moves with time. The fixed points are ( 13.171 ) for FPl-t, while 

(4 ±) , 2/2 ±} ) = H 1 , !) from E3S for FP2 ±- We also focus 

upon FP2±. 

A basis for the required analysis is prepared first by drawing 
the "null curves" for x' = and y' = 0. By choosing the 
vanishing left-hand sides of ( 13.511 ) and ( 13.521 ). we find that the 



former curves, solid (blue), are in fact hyperboloids, whereas 
the latter ones, dashed (red) curves are either hyperbolic or 
elliptic depending on £ 2 > 1/3 or C 2 < 1/3, respectively, as 
shown in Figs. [3]and|4j used separately for the analyses of the 
two choices for £ 2 . The crossings represent fixed points. Note 
that there are four fixed points (FP1±, FP2±) for C 2 < 1/3, 
while two fixed points (FP2±) for ( 2 > 1/3. 




FIG. 3: An example of £ 2 > 1/3 in the phase diagram of x = a' and y = In the overall diagram (a), we chose an example of 

£ 2 = 1/2 to draw solid (blue) and dashed (red) curves for the null curves for x' = and y = 0, respectively. They separate regions of 
different signs of them, denoted by + x and + y , for x' > and y > 0, respectively, for example. The crossings, or fixed points, occur between 
two hyperboloids, marked with a blob at FP2 + (a; = y = 1) as an attractor, and the one at FP2_(a; = y = — 1) as a repeller (no fixed points 
FP1± for this parameter). In each of smaller regions bounded by the null curves, as shown in (b) and (c), close-up views near the attractor and 
the repeller, respectively, we have directions of "flows" of solutions, or trajectories, shown symbolically by the arrows together with nearby 
pair of signs for the increment (+) or decrement (— ) of the x- and y-components denoted in this order. A dotted curve in (b) is a trajectory 
starting from the initial values lnp = 0.1, <p = 0.1, tp = at the initial Jordan frame "time" lnt = 1, where ip = (£/2)0 2 with £ = 1/4. 
The trajectory flows basically in accordance with the arrows, entering the frame of diagram near the lower-left corner, traversing to go outside 
beyond the right edge, re-entering at the top, spiraling finally into the attractor. 



In Fig. [3] (a), we draw two sets of hyperboloids for ( 2 > 
1/3. It is rather easy to determine which side of each curve 
implies positive or negative x' and y', as shown by the sym- 
bols like + x and + y , for x' > and y' > 0, respectively. 
We then determine in what direction a point, or better called 
a trajectory, "flows" with time in a given region in x-y space 
bounded by the null curves, as illustrated symbolically by ar- 
rows (green). In Fig. [3] (b), we show a typical trajectory 
represented by a dotted curve which enters the diagram first 
near the lower-left corner, making a big loop outside the frame 
of the diagram, then re-entering again, finally spirals into the 
crossing at FP2 + (x = y = 1), which corresponds to ( 13.341 ) 
and ( 13.351 ). This is the way we now establish our previous nu- 
merical results obtained on the basis of the heuristic approach 
G|], to be an authentic attractor in a strict sense. 

We notice, on the other hand, that a trajectory, or the so- 
lution, may not always converge to a fixed point, straying 
instead toward infinity, as will be illustrated in Appendix B. 
Obviously, however, not reaching the point of x' = y' = 
implies the destination not corresponding to the steady and 



lasting solution, as given by ( 13.341 ) and ( 13.351 ), for example. In 
other words, any solution that survives a long time must come 
from the attractor. 

There is another fixed point FP2_ (x = y = —1) to which 
no trajectory flows into as long as we start with a positive p, 
the same sign as Vq. Note that this fixed point corresponds 
to the contracting universe. The flows shown in Fig. [3] (c) 
indicates that this is a repeller to be interpreted as the time- 
reversed point against the attractor at FP2 + . 

We encounter another complication, on the other hand, if 
C 2 < 1/3, for which we have four crossings as illustrated 
in Fig. H|a). Focusing upon the behaviors in the upper-right 
quadrant, we have magnified views at each of the two, one in 
(b) around x = y = 1, corresponding to the fixed-point solu- 
tion of FP2 given by d3.34t -( [335] ) above, and another shown 
in (c) now categorized into the type of FP1. The flows in (b) 
do indicate that the the solution in all directions tends finally 
to the crossing, a real attractor, while according to those illus- 
trated in (c) there are some narrow strips sandwiched between 
two null curves in which the trajectory tends away from the 
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crossing. This reminds us of a saddle-point potential rather 
than the purely attractive or repulsive potential in mechanics, 
hence implying an unstable solution which may not survive a 
long time, leading ultimately to the same drift toward infinity, 
as was remarked toward the end of the discussion of the pre- 
ceding example in Fig. [3] For 1/4 < £ 2 < 1/3, therefore, 
we may fail to reach the attractor solution with some proba- 



bility depending on the initial values, or on what portion in 
the x-y plane we started off. We may constrain ourselves to 
C 2 > 1/3, though another a posteriori attitude might be sug- 
gested: Given what we are at present, right initial values, or 
right initial locations in phase space, must have been chosen 
to reach the attractor solution (|3.34| >-( f3735l >. 




FIG. 4: An example of 1/4 < ( 2 < 1/3. In the overall diagram (a), we chose an example of ( 2 = 0.2916, £ = 0.3890. Unlike in Fig. [3] 
the dashed (red) null curve is an ellipsoid, thus producing four crossings. In addition to the attractor at FP2+(:r = y = 1), denoted by a blob 
(green), we have another at FP1+ (x = 1.322, y — 1.134) marked by a cross (green), both accompanied with the time-reversal counterparts 
in the left-lower quadrant (FP1_ and FP2_), which we ignore for brevity. In a close-up view (b) around FP2+, we show an example of a 
trajectory with the same initial values as considered in Fig.[3](b), certainly spiraling into the fixed point in accordance with the arrows of flows, 
while we find no trajectories in (c) around another crossing of the type FP1 + . The behaviors of the arrows in the right-upper and left-lower 
strips remind us of a saddle-point potential resulting in unstable motion. 



We apply the same analysis for any values of 7, finding the 
similar results. We summarize the power index of scale factor 
of the the attractor solutions (FP1 + for £ < 1/2, and FP2 + 
for ( > 1/2) in Fig. [5] We find that the power-law inflation 
appears when Q 2 < 1/8. The power exponent of the scale 
factor in the Einstein frame is always 1/2 for FP2 + , which 
does not depend on the equation of state of the matter fluid, as 
derived in (13.35b . 

IV. COSMOLOGY IN THE JORDAN FRAME 

We repeat the similar analysis for the same system as in 
Section III now in the Jordan frame. 

From (12. U together with the FRW metric 

ds 2 = -dt 2 + a 2 dsj , (4.1) 

we derive 

3£0 2 (h 2 + + 6£H<M> =^<p 2 + V a +p (4.2) 

e(00 + 2 + 3w) = C 2 (4^o + P - 3P) (4.3) 
p + 3H(P + p) = (4.4) 



where 

H=-, (4.5) 
a 

while P = P* cxp(4£er) and p = p* exp(4£er) are the pres- 
sure and the energy density, respectively, in the Jordan frame. 
Note that the quantities in the Jordan frame are all denoted 
unmarked in contrast to those marked with the subscript * in 
the Einstein frame. Accordingly, the dots in fl4.2| )-( |4~4l imply 
differentiation with respect to t instead of t*. 

Since ds 2 = Vt~ 2 ds 2 with Q 2 = £0 2 = exp(2Cc), we have 
the relations between the variables in Jordan frame and those 
in the Einstein frame as 

dt = n~ 1 dU, (4.6) 
a = fT 1 ^ . (4.7) 

Note that we do not change the coordinate system when we 
perform a conformal transformation. However, since we use a 
cosmic time in each frame, we have to change the time coordi- 
nate between i* and t according to d4.61 l. The Hubble expan- 
sion parameters H and iJ* are defined by each cosmic time 
as 

a at a* ctr* 
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FIG. 5: The power index p* of the scale factor in the Einstein frame 
for the attractor solutions. The power-law inflation appears when 
< 1/8. The power exponent is always 1/2 for FP2 + , independent 
of 7. 



Hence we have the relation 



This is again a self-autonomous system with two variables x 
and y: 



3C 
4 
1 



[2-fx 2 -4xy + 4( 2 (4-3j)y 2 + 7] , (4.18) 



y = — [-2(2 - 7)^ + \% l xy - 12 7 CV + 7 • 

(4.19) 

where x = C,^ 1 ®' and y = J£ . These equations turn out to 
be precisely the same as Eqs. d3.15l l and ( 13.161 ). respectively, 
implying the same dynamical system, sharing the same fixed 
points: FP1± and FP2+ d3~T8l 
The energy density is given by 



P_ 
To 



12IK 2 



r(*') 



' \2 



= l2Q 2 y 2 - 2x 2 - 1 , 



(4.20) 



precisely the same as (13.191 l for the Einstein frame. 

We only show the explicit solutions of the fixed points, be- 
cause the dynamical properties such as an attractor is the same 
as that in the Einstein frame. For FP1±, we have 



H * — fi 



H 



din ft 

dt 



(4.9) 



Introducing new time coordinate r and and new scalar field 
$, which are defined by 



(It = 2C 1/2 V Q 1/2 cj)- 1 dt 
<E> = In (A, 



respectively, we rewrite Eqs. ( 14.2l )-( l4~4"b as 

3/c£e 2<E> 1 „„, 1 



63i 2 



— ($') 2 T 

2V a? C 2 2 



where 



$" _ 2($') 2 + 'iK®' = C 2 



a 



(4 - 3 7 )p 

4Vn 



(4.10) 
(4.11) 



(4.12) 

(4.13) 
(4.14) 



(4.15) 



_1 

a = a (T-T ) 4 C 
2C 



2-1 



±- 



(r-ro) 



(4.21) 



p = 0. 



This gives a power-law inflation for £ 2 < 1/ 
For FP2± , we have 



±- 



V7 



P = 



^2(2- 7 -2(4-37)0 
2(4C 2 - 1) 



2 - 7 - 2(4 - 3 7 )C : 



V , 



(4.22) 
(r-T ), (4.23) 

(4.24) 



where ao and To are integration constants. The space time is a 
static Minkowski space, in contrast with the expansion in the 
Einstein frame as shown by ( 13.351 1. 

The vacuum solution in radiation-dominance in ^ may be 
interpreted as the limit C 2 — > 1/4 in (I4.22b - d4.24l l. The con- 
stant a was also derived in lfl8ll . 



and the prime is the derivative with respect to r. We also 
assume the equation of state, P = (7 — I) p. 

Now, we consider only the spatially flat universe, i.e, k = 
(See Appendix A for k 7^ 0). Taking the derivative of Eq. 
(l4~T2l and using Eqs. (l4~T2l , gjj) and (l4~14l) . we find 



IK' = ^tt^ ($')' + 2^$' - ylK 2 



4C 2 

2 



(4.16) 



$" = ^($') 2 - 3W$' + 3(4 - 3 7 )C 2 ^ 2 + ^C 2 



(4.17) 



V. COSMOLOGY WITH A POWER-LAW POTENTIAL 

The analysis in the preceding section can be readily ex- 
tended to include the power-law potential, 



v(<f>) = x a <p a 



(5.1) 



in the Jordan frame. An example is provided by extending the 
action ( |2.5t in string theory to 



-2<I> 



R(g) - 2A + 4(V$)' 



, (5.2) 
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in which we have added the term of A, to be corresponded to 
d5.lt by choosing a = 2 and A 2 = A/4. 

The action equivalent to ( 15. Q but expressed in the Einstein 
frame is given by 



S = J d 4 Xy^j 

+ / d 4 x^/-g*L(a, ip*,g*) 



ii?,(^)-I(V*a) 2 -V» 



where 



V* = exp[-C(4 - a)a] V , 



(53) 



(5.4) 



The basic equations for cosmology are then; 
k 1 

#* + ^ = - (<x 2 + K +/>*), 
dV* 

a + 3H,& + = ((p* - 3P*), 
do 

p* + 3-ff*(P* + p*) = -(<r(p* — 3P„). 
Introducing a new dimension-free time coordinate as 

dr* = 2\/V«dt.. 1 



(5.5) 

(5.6) 
(5.7) 



(5.8) 



with Vq = A Q /|£| a / 2 . Since this is of the same form as Eq. 
( 12.91 . we repeat the same analysis as before. One of the points but replace (| 3.13 |) by 
to be kept in mind is that the conformal transformation is al- 
ways the same as (12.6b with (12.7b . Hence we have the same 
energy-momentum conservation law (13.41 1. as well. 



also assuming the equation of state and focusing on k = as 
before, we obtain formally the same equations as (13.91 l-( l3~TTl ). 



= - [3 7 Jf* - C(3 7 - a)a'] . (5.9) 



Computing in the same way as before, we come to replacing ( 13.151 ) and ( 13.161 ) by 



x' = | [2 (37 - a) x 2 + 12C 2 (4 - 3-f)y 2 - 12xy + 3 7 - a] 
y' = ^ [-2(2 - !)x 2 - 12 7 CV + 4C 2 (4 - a) xy + 7 ] , 



respectively, where x = a' and y = ( 1 J{* as before. 
We find four fixed points as the previous case: 



FPl ± :(x F ,y F ) = (x[ ± \y[ ±) )=± 



FP2± : (x F ,y F ) = (x 2 ±) ,y { 2 ±) )=± 



1 ( C(4 - a), - 

v/2(6-C 2 (4-a) 2 ) V C 

1 



v/6[3 7 (2- 7 )-2(4-3 7 )(3 7 -a)C 2 ; 
I 



(3 7 , 3 7 - a) 



(5.10) 
(5.11) 



(5.12) 
(5.13) 



The fixed points FP1± exist if 



(5.14) 



For fixed points FP2±, the constraint for existence is a little 
more complicated, as we find; 



a > 3 7 



1 - 



2- 7 1 



2(4 - 3 7 ) C 2 
any values of £ and a 

2-7 1 



a < 37 



1 



2(3 7 - 4) C 2 



for 7 < - (5.15) 
for 7 = \ (5.16) 
for 7> \- (5.17) 



Next we analyze the solutions of these fixed points. In what 
follows, we consider only the expanding universe (y F > 0), 
i.e. FP1 + , while FP2 + for a < 3j and FP2_ for a > 37. 



For FP1 + , we find 



C(4 - a) 



^2(6-C 2 (4 -a) 2 ) 

T* 



T* + CO 



a* = a^o exp 



V2(6-C 2 (4-a) 2 ) 
The cosmic time t» is given by 

[ C 2 (4-a) 2 n 
i* = £*o cxp — . 

P [^2(6-C 2 (4-«) 2 ) 

Hence the solution of FP1+ is described by 



2 . ( U 
In 



C(4-a) V** 



CO 



(5.18) 
(5.19) 



(5.20) 



(5.21) 
(5.22) 
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where 



P* 



C 2 (4-a) 2 



(5.23) 



In this case, matter does not contribute to the expansion of the 
universe, i.e. at the fixed point, = 0. The cosmic expansion 
becomes inflationary if 



(5.24) 



For FP2±, we have 
a = ±- 



37 T* 



v /6[3 7 (2-7)-2(4-3 7 )(37-a)C 2 ] 
±(37 - a)Cr* 



(5.25) 



a*o exp 



x /6[3 7 (2-7)-2(4-3 7 )(3 7 -a)C 2 ] 

(5.26) 



The cosmic time i* is given by 

±3 7 (4C - a) r* 



£* = i*o exp 



2 V6[3 7 (2 - 7) - 2(4 - 3 7 )(3 7 - a)C 2 



(5.27) 



Hence the solution of FP2+ is described by 



In 



C(4-a) V**o 
f 



where 



where 



2(3 7 -a) 2 
2(3 7 -4) 



Ap* = 



37(4- a) ' 



(5.28) 
(5.29) 

(5.30) 

(5.31) 



which describes the deviation from conventional power expo- 
nent with the adiabatic index 7. It is due to the interaction 
between the matter fluid and the scalar field (f>. It is precisely 
this interaction that keeps the ratio p*/V* constant at FP2±. 
The value is obtained from Eq. ( 13. 19b ; 



6 W* 



2C 2 « ±: 



FP2 



2[(4-a) 2 + (37-4)(4-a)-3 7 /C 2 ] 
2(3 7 - 4) 2 + 2(3 7 - 4) (4 - a) + 3 7 (2 - 7 )/C 2 ' 

(5.32) 



which is consistent with d5.9t for the energy density. 

In order for this energy density to be positive, we have to 
impose the following condition: 



1 

a < — 
2 



(3 7 + 4) - v/(3 7 - 4) 2 + 12 7 /C 2 for FP2+ 



a > 



(3 7 + 4) + ^(37 - 4) 2 + I27/C 2 



for FP2_ 



We also find a simple result = 1/2 follows either by 
a = for any 7 or by 7 = 4/3 for any a. This result may be 
described by 



where 



Apt 



P* = - + Aj/„ 



(37 - 4) a = a 
67(4 -a) 4 P * 



(5.33) 



(5.34) 



We are now looking into more details of the power-law in- 
flation. We find a power-law inflation by a if the power expo- 
nent of the potential V* is sufficiently small, i.e. |£(4 — a)| < 
V2, and if there is no coupling between the matter fluid and 
the scalar field a in the Einstein frame. This type of inflation 
is realized in the fixed point FP1 + even when we include the 
coupling with the matter fluid, which has nothing to do with 
inflation at FP1 + . However, the power-law inflation may oc- 
cur also for FP2±, in which the coupling to the matter fluid 
is non-negligible. From the condition > 1 for the power 
exponent of the scale factor in the Einstein frame, we find the 
following conditions: 



a < 



37-2 



a > 4 



37-2 



for 


7< i 


(5.35) 


for 


2 

7= 3 


(5.36) 


for 


2 

3<7< 


^(5.37) 


for 


4 

7= 3 


(5.38) 


for 


4 

7> 3- 


(5.39) 



4 < a < 

no case 

67 
37-2 ' 

This is a new type of inflation. The potential itself is too 
steep to cause inflation, but the matter fluid assist to cause 
a faster expansion because of its coupling to the scalar field. 
Note that including the matter fluid with 7 < 2/3 means that 
we can have inflation just by this fluid. But the coupling with 
the scalar field does not assist such a type of inflation. Rather 
it will restrict the possibility of inflationary expansion. 

Next we have to analyze the stability of the fixed points. 
We can analyze it by perturbations or by use of the phase di- 
agram, as shown in the previous section. Here we give the 
perturbation analysis. Inserting the perturbations around the 
fixed points x = Xf + Sx, y = ijf + Sy in the basic equa- 
tions ( I5.10l l and ( 15.111 ). we find a set of the linear perturbation 



11 



equations ( 13.391 ) with the components of the matrix being (13.41b with 

' A xx = (37 - a) x F - 3y F 
A xy = -3x F + 6C 2 (4 - 3 7 )y F 

A 

A 



2-7 . (4-q) 



(5.40) 



(4-a) , 
>2 XF ~ 3"fVF 



Setting Sx,Sy oc e UT * , we find the eigen equation for u> as 



J 



TiA = (3 7 + 2-3a/2)a; F -3(7 + l)yF (5.41) 
det.4 - ^ - a)(4 - a)( 2 - 3(2 ~ a)] x 2 F 



+ 3(aj — IO7 + 8)x F y F 

+ 3[3 7 -C 2 (4-3 7 )(4-a)]?,|. 



(5.42) 




FP2-I 





FP2-I 



(a) 7 = 1 (dust) 



(b) 7 = 4/3 (radiation) 



(c) 7 = 2 (stiff matter) 



FIG. 6: For various matter fluids ((a) dust [7 = 1], (b) radiation [7 = 4/3], (c) stiff matter [7 = 2]), we show the ranges in £-a plane 
of attractor fixed points (FP1 + ,FP2±). The shaded regions FP1 and FP1-I are those for FP1 + , while FP2 and FP2-I are those for FP2±, 
respectively. The dark shaded regions FP1-I and FP2-I are the ranges for power-law inflation.There is no attractor solution in the unshaded 
region. 



Using these equations, we analyze the stability of the fixed 
points as before. The result is shown in Fig. [6] The regions 
denoted by FP1 and FP1-I give the attractor solution of FP1+, 
while those by FP2 and FP2-I correspond to the attractors of 
FP2±. ForFPl + , the inflationary solutions exist in the case of 
a flat potential (C|4 — a\ < 1/2), but the new type of inflation- 
ary solution also appears for FP2± if 7 7^ 4/3. Each region 
for inflation is also shown by FP1-I or FP2-I. This new type 
of inflation is assisted by a coupling between the scalar field 
a and the matter fluid. 



In the Jordan frame, we show only the behavior of the fixed 
points FP1± and FP2±. After a conformal transformation, we 
find them: 



where 



a oc t p 
4> oc t q . 



(5.43) 
(5.44) 



2(1 - (4 -q)C 2 ) 2 

(2-a)(4-a)C 2 '^ 

2a 2 \ 
2 - a' 2 - a J 



for FP1± 

for FP2± 

(5.45) 



The scale factor in the Jordan frame for FP2± is not con- 
stant except for a = 0. It is expanding for a < or a > 2, 
while it is contracting for < a < 2. For a = 2, we find 

a oc exp[— 



VI. CONCLUDING REMARKS 

We have presented a formulation in which we trace the 
temporal development of cosmological solutions of the scalar- 
tensor theory in two-dimensional phase space. Thanks to as- 
suming a simple equation of state, in Sections III, IV and V, 
we have obtained two different sets of fixed points, FP1± and 
FP2±. Conditions of obtaining attractors are studied in detail. 
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We have established the attractor nature of the fixed points. 

1 12 

At FP2±, the scale factor behaving like a constant and tj in 
the Jordan and the Einstein frames, respectively, for C, 2 > 1/4 
or £ < 1/2 (u> < —1/2) with e = — 1, when we have a simple 
cosmological constant in the Jordan frame. This solution is 
also accompanied with the proportionality between and V*, 
called a scaling behavior, which is going to be replaced by the 
"interlacing" behavior, as exemplified in Fig. 5.8. of y|], by 
further extending the model. 

An extension to the power-law potential in Section V has 
shown that the coupling between the scalar field and the matter 
provides a new type of inflation with 7 off the conventional 
choice 4/3, even if the potential of the scalar field is too steep 
to cause inflationary expansion by itself. 

We have also learned that, even confining ourselves to the 
fixed-point solutions, reaching an attractor can be somewhat 
complicated if we have another fixed point, as we faced in 
the example of Fig. 21 though a simple recipe is shown to 
be applied by adopting either of the two attitudes suggested; 
limiting the range of £ 2 or being content with an a posteriori 
consideration that initial values of x and y had happened to be 
in favor of reaching the attractor. The first appears to be the 
case in our realistic choice (( = 1.5823) to fit the observed ac- 
celerating universe shown in Fig. 5.8 of fH, for example. The 
second of the above suggestions might also apply to drifting 
to infinities, as will be demonstrated in Appendix B. 

If there is a spatial curvature, the result will be changed. 
However, we can show that the fixed point is still an attractor 
if the universe expands very fast, i.e, if it is an inflationary 
universe, as shown in Appendix A. 

The most intriguing result in the present scalar-tensor the- 
ory is, however, that the static universe in the Jordan frame 
is an attractor, an unavoidable fate in the presence of the cos- 



mological constant. There is no smooth limit as Vq — > 0. As 
we re-iterate, the Jordan frame features truly constant masses 
of microscopic fields, according to the Brans-Dicke model, 
originally intended to qualify this frame to be physical, allow- 
ing a non-static universe in the absence of the cosmological 
constant. Its presence alters the entire situation, forcing us to 
accept ma = const. This entails eventually that the universe 
in the Einstein frame expands in the same rate as the (time- 
dependent) microscopic length standard. This crisis will be 
evaded only by leaving the Brans-Dicke model, as was elabo- 
rated in Jltljlllll together with the ensuing consequences. The 
view that this crisis hinges upon the attractor nature of the so- 
lutions is now reinforced even more strongly by our study in 
this article. 

We also emphasize that the extension to the power-law po- 
tential leaves the above crisis unsolved. The exponent 1/2 in 
radiation-dominance in the Einstein frame is a unique conse- 
quence of 7 = 4/3 independent of the way Vq is modified by 
the scalar field. The argument on the power-law inflation also 
remains unaffected, as we point out, by the structure of the 
mass term. 
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APPENDIX A: THE EFFECT OF CURVATURE TERM 

(MO) 

In this Appendix, we study the curvature effect. It may be 
convenient to analyze the equations rewritten by new variables 
because the fixed points are constant. We discuss the cases 
with a cosmological constant and with a power-law potential 
separately. 



1. The case with a cosmological constant 

We show the curvature effects both in the Einstein and in 
the Jordan frames in this order. 



a. Curvature term in the Einstein frame 

The curvature term is proportional to e^/a 2 . We evaluate 
the time evolution of this term by 



2. The case with a power-law potential 

In this case, we can repeat the same analysis. The curvature 
term is proportional to e' 4_Q ^ CT /a 2 . So the time evolution 
near the fixed point is given by 



dln(e 



(4-a)C«r/ 2) 



= C [2{2xf - Vf) - ccxf] 



2C 2 (2 - a/2) 2 - 1 

CV2(6-C 2 («-4C) 2 

C[6 7 -a(37-2)] 

V6[3 7 (2-7)- 2(4- 3 7 )(37-a)< 2 ] 



for FP1+ 

for FP2+ 
(A5) 



The curvature term is not important for inflationary solution 
near the fixed point FP1 + (C 2 (4 — a) 2 < 2) and near the 
fixed point FP2+ (a > 67/(37 - 2) with 7 > 4/3) or FP2_ 
(a < 67/(37 - 2) with 2/3 < 7 < 4/3). 



dr* 



= 2{2x - y) 



(Al) 



We find the behavior near the fixed point (xp, Vf) as 
d\n{e<°/al) 



dr* 



= 2{2x F - Vf) 



= < 



8C 2 



2<V3 - 8C 2 



for FP1. 



V7 



(A2) 



{ V2(2- 7 -2(4-3 7 )C 2 ) 



for FP2 



+ 



Hence we find that the curvature term decreases in time for 
inflationary solution at FP1 + (£ 2 < 1/8). It can be ignored. 
However, it will grow in time if the universe expands without 
acceleration. It will becomes important as the same as the 
usual case. 



b. Curvature term in the Jordan frame 

The time evolution of the curvature term in the Jordan 
frame is given by 



d\u{4> 2 /a 2 ) 



= 2($' - <j>M) = 2(2$' - JC) = 2{2x - y) 



We find the behavior near the fixed point (x f,Vf) as 
rfln(0 2 /a 2 ) 



dr 



= 2{2x F - y F ) , 



(A3) 



(A4) 



which is exactly the same as Eq. ( IA2I ). Hence the curvature 
term near the fixed point FP1 + is not important for inflation- 
ary solution (C 2 < 1/8). 



APPENDIX B: ADDITIONAL SOLUTIONS FOR 
RADIATION-DOMINANCE IN THE JORDAN FRAME 

We often relied on the numerical approach to solve the cos- 
mological equations in the Jordan frame because the solution 
is characterized by the simplest aspect of the static universe. 
We encounter, however, another complicated aspect to be dis- 
cussed in what follows. 

We started conveniently from 



6ipH 2 
p + AHp 



2 V 

0, 



Vo + p- QHlp, 



(Bl) 

(B2) 
(B3) 



where tp = (£/2)</> 2 in terms of which (14.3b has been put into a 
simplified form in dB21 >, as in J^. We also write V$ = A = 1. 

Since H occurs always without derivative, we may elimi- 
nate it by using (IB3t . for example; 



H = (B4) 

4p 

Eqs. ( IBU and (IB2t are then put into 

/ • \ 2 • -2 

3(£) ^-12^ + 2r J — =8 (Vo + p), (B5) 



Ap 



(B6) 



which are to be solved by giving three initial values of p, tp, tp. 

We notice, however, that we solve dB5b and dB6b first 
with respect to p/p. This involves solving an algebraically 
quadratic equation for p/p, thus producing two differential 
equations, hence resulting in two separate solutions. An ex- 
ample of numerical solutions is shown in Figs. |7fa)and(6). In 
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Fig. |7](a), developed basically from Fig. 4. 1 of (3f|, we do find 
asymptotic behaviors H — > 0, 4> — > \/4Vb/(6£ — 1), p — > 
— 3Vo(2£ — l)/(6£ — 1) corresponding to an attractor solu- 
tion. Fig. PTJb) illustrates, on the other hand, the solution of 
another equation, but sharing the same initial values of p, ip 



and ip as discussed in (a). This one represents, however, a 
shrinking universe taking place in a short time. This type of 
the second solution occurs nearly always. It even appears as 
if we are going to lose an opportunity to reach the fixed-point 
attractor. 




5 10 15 20 0.05 0.1 0.15 0.2 1 2 3 4 



In t In t x 

FIG. 7: Two solutions for the initial values at lnt = 0; p = 0.1, tp — 0.1, <p — together with Vo = 1. The asymptotic values mentioned 
above, hence d3.34t -l P3.35t are reached in (a), while these quantities in (b) diverge quickly at lnt — + 0.248, in this example. Behaviors of the 
solutions shown in (a) and (b) are represented by two different trajectories in phase space of x, y of (c), starting at two different points A and 
B. The upper curve corresponding to (a) spirals finally down into the attractor denoted by a blob (green) at x — y — 1, while the lower curve 
for (b) runs steadily toward the right-lower infinity, as "guided" by the arrow. 

I 



Fortunately, as it turns out however, the phase-space de- 
scription is so generous that the resulting two solutions are 
accommodated as those due to two different initial locations 
of y, as demonstrated in Fig. |3 C )- Again for a typical solution 
for £ = 1/4, we find two trajectories for the solutions of < IB5b - 
dB61 l. Both start at x = corresponding to our initial value 
ip' — 0, but with y separated from y — nearly the same dis- 
tance in the opposite directions. The upper curve starts at the 
point A. After a loop-like trip, it finally spirals down into the 
fixed point at x = y = 1, as usual. Another curve starting at 
B drifts steadily toward x — ► oo, y — ► — oo. Both trajectories 
started immediately "outside" the y' = boundaries. 

As we find from y = preceding d3.15t and (13 . 1 6b . 

the solutions with y > and y < correspond to the ex- 
panding and contracting universes, respectively, in the Ein- 
stein frame. From this point of view, what we are in the ex- 
panding universe at quite a late epoch must be a consequence 
of the initial value selected to be y > 0, as represented by the 
upper trajectory in Fig. |7](c). In this sense the solution like 
Fig. |7](b) is excluded leaving the attractor solution alone. 



APPENDIX C: DUST-DOMINANCE IN THE 
SCALE-INVARIANT MODEL 

In sections III and IV we discussed cosmological solu- 
tions mainly in the radiation-dominated universe finding a 
crisis arising from too much time-dependent masses of parti- 
cles evaded finally by departing from the Brans-Dicke model, 



even at the risk of WEP violation. The same type of analy- 
sis of dust-dominance suffers more seriously because it en- 
tails a* ~ t x J 2 , as shown by d3.351 l even for 7=1. As 
was discussed in Oft, the remedy comes simultaneously from 
the scale-invariant model intended to overcome the crisis for 
radiation-dominance. We sketch below how this model pro- 
vides attractor solutions also for the dust-dominated universe. 
See Chapter 4.4.3 of |H and Section 3.4 of H for more de- 
tails. 

The field equations for a and p turn out to be given by d3.6t 
and fl3.7b with the right-hand sides removed to the classical 
approximation. 

As a remarkable difference from d3.35t we find 

a* = a»o (f* A*o) 2 ^ 3 1 (CI) 

in agreement with the conventional law of expansion. Eqs. 
( 13.51| l and d3.52t are replaced by 

x' = C (2a: 2 - 3xy + l) , (C2) 

V' = ^ [~\x 2 + \ + ( 2 (ix - 3y)y^j , (C3) 

respectively. The solutions with (13.341 ). ( IC1I ) and 

exp(-4c> ) = i;C 2 , (C4) 
16 

in place of d3.37| ) are obtained for x = l/\/2, y = 2\/2/3, 
corresponding to an attractor yielding x' = y' = 0. 
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A similar distinction between the elliptic and the hyperbolic 
curves as in radiation-dominance occurs also for £ 2 < 3/8 
and ( 2 > 3/8, respectively. The same recipe should apply as 
mentioned toward the end of section VI. 

We add that the scale invariance coming from the absence 
of dimensional coupling constants has an advantage that a 
serves as a massless Nambu-Goldstone boson which will ac- 
quire a small mass after the invariance is finally broken ex- 
plicitly through loops, as discussed in Section 6.3 of 



APPENDIX D: ANOTHER APPROACH TO THE 
POWER-LAW POTENTIAL 

It seems also useful to apply (IB 1 b - (lB3b to the power-law 
potential to offer a simplified alternative to derive the same 
common result = 1/2 for radiation-dominance as stated 
immediately after in (15.301 . 

We multiply V in (IBU and ( |B3t by (j) a . We search for the 
solution of the type 



and 



»(*) ~ t P , 



4>{t) ~ 



(Dl) 



(D2) 



By substituting them into (IB2t modified as above and com- 
paring the exponents of t we obtain 







(D3) 



implying that a = corresponds to (3 = 1. In the modified 
( IB 11 1, on the other hand, we find all the terms other than p 
to behave like t 2 °~ 2 , while dB3b entails p ~ t~ Ap . For a 
consistent approach we expect 2/3 — 2 = — 4p, or 



2(2- a)' 



(D4) 



This point was not properly recognized when it was erro- 
neously stated in Appendix B of j5] that only a = is consis- 
tent with the physically acceptable condition p* = 1/2 (The 
coefficients a and j3 in Appendix B of ^ are replaced by a/ 2 
and (3/2, respectively, according to the present notation.). 



Now from (14.6b and Q ~ <fs combined with ( ID2b . we obtain 

dt* = Qdt ~ tPdt, (D5) 



which is integrated to give 



(D6) 



where we have ignored inessential coefficients for simplicity. 
In the same context we also use d4.7l i to derive 



(D7) 



Combining this with ( ID4l i and ( ID6t we identify the right-hand 
side with tZ* where 



p* = 



(D8) 



which turns out to be the same as the result for the purely 
constant Vq = A, in agreement with ( 15.301 with a = for any 
7. This justifies that the present solution to be an attractor. 

We also notice that the universe is no longer static in the 
Jordan frame, as shown in ( ID4I ). We may no longer rely on 
the simplest argument am = a^m* = const to leave the BD 
model. According to Appendix D of [5], particularly its (3.10) 
and a more general procedure developed there, however, the 
mass m* of the matter fields in the Einstein frame is related to 
m in the Jordan frame as 

= fT 1 ™ - t~ m. (D9) 

From ( ID3b and dD6l ) we find 

t -P „ t -P/(ft+l) „ c (l/2)/(l-a/4) _ (D10) 

The exponent — (l/2)/(l — a/4) is not exactly the same as 
— 1/2 which would have implied that the universe looks static 
if measured with respect to the microscopic length standard, 
but is nevertheless far from zero as expected if the Einstein 
frame is qualified to be a physical conformal frame for any 
reasonable choice of a. In this sense the crisis for the purely 
constant Vq as discussed before is not evaded by multiplying 
it by the scalar field. Departure from the BD model seems still 
unavoidable. We also add that the argument for the constant 
ra* with the assumed scale-invariant model remains unaltered 
by the multiplied scalar field. 



